Two-Dimensional Quantum XY Model with Ring Exchange and External Field 
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We present the zero-temperature phase diagram of a square lattice quantum spin 1/2 XY model 
with four-site ring exchange in a uniform external magnetic field. Using quantum Monte Carlo 
techniques, we identify various quantum phase transitions between the XY-order, striped or valence 
bond solid, staggered Neel antiferromagnet and fully polarized ground states of the model. We find 
no evidence for a quantum spin liquid phase. 



Studies of two-dimensional spin-1/2 quantum magnet 
and boson models have provided insight into novel quan- 
tum phases and quantum critical points Q. Recently, 
interest has focused on models which have multi-site ring 
exchange The ring exchange interaction, either 

alone or in competition with the usual spin or boson near- 
neighbor exchange, has been shown to promote a variety 
of exotic quantum ground states Q, including in some 
cases a spin-liquid state 0. Of particular importance 
is the class of two-dimensional model Hamiltonians that 
contain quantum spin-1/2 or boson operators interacting 
with ring exchange that can be simulated using quan- 
tum Monte Carlo (QMC) techniques without a negative 
sign problem. With modern algorithms, such models can 
be studied numerically on large lattices without approx- 
imation, providing a laboratory for surveying the critical 
behavior that separates various quantum phases. 

One important model in this respect is the easy-plane 
J-K model Js] that has quantum S = 1/2 spins on a 
square lattice with a near-neighbor exchange J and a 
four-site ring exchange K. This Hamiltonian is par- 
tially motivated by the undoped cuprate materials [^, 
where ring-exchange processes are believed to contribute 
to experimental signatures beyond those explained by the 
near-neighbor Heisenberg model. The two-parameter J- 
K model, despite its simplicity, displays a surprisingly 
rich and complex phase diagram (sl , with three dis- 
tinct zero-temperature phases. These are an XY-ordered 
or superfluid phase for large J, a staggered Neel or bo- 
son charge density wave (CDW) phase for large K, and 
a striped or valence bond solid (VBS) phase for interme- 
diate K/ J . The zero-temperature phase transition be- 
tween the VBS and Neel phases is first order, however 
previous numerical results 5, 7] indicate the existence of 
a continuous quantum critical point (QCP) at the zero- 
temperature superfluid- VBS boundary. 

The question naturally arises as to the behavior of the 
easy-plane J-K model under the influence of a magnetic 
field. This is interesting both as a study of the evolution 
of the QCP, as well as the behavior of the ground state 
phases away from half-filling. Using stochastic series ex- 
pansion (SSE) QMC, we present here the basic features 
of the zcro-tcmperature phase diagram of the easy-plane 
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FIG. 1; The schematic zero-temperature phase diagram of 
the easy-plane J-K-h model. Phase boundaries are drawn as 
solid lines. Dashed lines indicate cuts along which we have 
examined the transitions between the various phases, as dis- 
cussed in the text. 



model described by the Hamiltonian 
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where Sf is the z component of a quantum spin-1/2, 
Bij = S^Sj' + Sj' is a near-neighbor exchange, and 
Pijki — SJ S'^ Sf + Sf Sj' Sf generates a four-site 
ring exchange. Here, (ij) denotes a pair of nearest- 
neighbor sites and (ijkl) are sites on the corners of a 
square plaquette on the L x L lattice. For K = 0, this is 
the standard XY-model in a uniform magnetic field, or 
alternatively hard-core bosons with a chemical potential. 
For h = 0, this model is in a XY-ordered or superfluid 
phase for temperatures less than the Kosterlitz-Thouless 
transition temperature of Tkt/J ~ 0.68 0- With the 
application of a uniform magnetic field the average mag- 
netization m ~ (S^) of the XY-superfluid increases from 
zero (m = 0) until it saturates into a fully spin polarized 
(to = 1/2) state at h/J = 4 For J = and ft, = 0, 
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the ground state of the system is Neel antiferromagnetic 
order 0, For /i = 0, it was found 0, that an in- 
termediate VBS phase exists for 7.9 < K/ J < 14.5, in 
which the expectation value {Piju) alternates in strength 
with a period of two lattice spacings in one of the lattice 
directions, suggesting the term "striped" order. 

To study the effect of the uniform magnetic field h 
on the ground state properties of the easy-plane J-K 
model, we use the SSE quantum Monte Carlo simulation 
method [lO| that was previously applied to the /i = 0, J- 
K modelja, • In order to implement the SSE method, 
the operators in the Hamiltonian are represented as 
four-spin plaquette operators. Diagonal operators involv- 
ing h terms are added to or removed from the SSE ba- 
sis state expansion using a simple Metropolis probabil- 
ity algorithm. Off-diagonal [J oi K term) operators are 
sampled using the directed-loop algorithm |3)|lOj, which 
becomes increasingly important for simulation efficiency 
with increasing magnetic field strength. The directed- 
loop equations for the J-K-h model are only slightly 
more complicated than for the pure J-K model 0, and 
are presented elsewhere flU . The QMC algorithms were 
tested on L = 4 lattice sizes against exact diagonalization 
results and previous QMC simulations on the pure XY 
and J-K models. In this paper, simulations were carried 
out on square lattices of linear dimension L (number of 
spins N — L^), at temperatures T — 1/ (3 low enough to 
ensure convergence into the ground state. 

A variety of physical observables of direct relevance 
to the ground states of the model are accessible through 
the SSE method. It is straightforward to calculate the 
internal energy [loj since its statistical estimator is just 
the number n of plaquette operators in the SSE basis- 
expansion operator sequence multiplied by T: E — 
— (n) I f]. The spin stiffness (or superfluid density in the 
boson representation) is defined in terms of the energy 
response to a twist cj) in the periodic boundary of the 
lattice by 



(2) 



and is directly estimated using the winding number fluc- 
tuations in the SSE simulation 12]. In addition we cal- 
culate the plaquette structure factor 
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Here, 
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ai, . . . , 04 are the sites belonging to plaquette a, 
located at r^. In the VBS phase, the square of the 
magnitude of the order parameter per site is {Mp) — 
[^^(Tr, 0) -I- S'p(0, 7r)]/2L^. Similarly, the square of the or- 
der parameter {Ms) of the Neel ordered phase is obtained 
from the structure factor 



ffc)-q 



(4) 



with (Ms)'^ = Ss{tt, 7r)/L^. Here, j and k are lattice sites 
located at lattice coordinate r^. The quantities {Mp)^ 
and (Ms) are expected to decrease as 1/L^ (signifying 
short-range correlations) in phases without the respective 
order, but tend to a finite value for large L in phases 
where long-range order occurs. 
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FIG. 2: The ground-state energy [E) and magnetization (m) 
of an L = 8 system along cut (A) in the phase diagram Fig.0 
This set of simulations was performed with parameters K/ J = 
14 and /3J = 320. The hysteresis effects were obtained by 
systematically increasing and then decreasing K/J in steps, 
with system configurations stored at the end of one K/J step 
and used to begin the next. 



By directly observing the behavior of the spin stiffness 
(superfluid density) and the VBS and Neel order param- 
eters, we are able to map out the phase boundaries of 
the J-K-h model as illustrated in Fig. ^ In general, we 
find no persistent regions of quantum disorder (i.e. a spin 
liquid state) in the vicinity of the h = quantum criti- 
cal point. Rather, the QCP appears to evolve smoothly 
into a quantum phase transition between the superfluid 
and VBS regions for < /i ^ 6. The J-K-h model also 
exhibits a direct superfluid to Neel order transition for 
6 ^ ^ 11, a feature not contained in the h = phase 
diagram. Finally, for large h, the model finds a fully po- 
larized spin state with m = 1/2. This latter phase tran- 
sition is strongly first order for K/J ^ 5, displaying pro- 
nounced metastability and hysteresis effects in the sim- 
ulation (see Fig. 121). Renormalization Group treatments 
of two-dimensional bosons [l^ , as well as spin- wave cor- 
rected mean-field theory and simulations of a hard-core 
boson Hamiltonian indicate that at K = 0, the pure 
XY model exhibits a continuous transition to the fully 
polarized state at h — 4J. This suggests that a tricriti- 
cal point (TCP in Fig. exists on the phase boundary 
somewhere between < K/J < 5, above which the tran- 
sition to the fully polarized state becomes first order. 
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The energy crossover and density hysteresis of Fig. [21 
provide one indicator of a first-order transition. Aher- 
natively, one may look for an abrupt discontinuity in the 
order parameter (for large system sizes) or for double- 
peaked probability histograms for data in the transition 
region. To illustrate this we turn now to a detailed set of 
simulation results for the superfluid-Neel phase bound- 
ary along cut (B) in Fig. Q As illustrated in Fig. |21 the 
boson and superfluid densities develop significant discon- 
tinuities for larger systems as the phase boundary is tra- 
versed. This abrupt discontinuity does not appear for 
L < 20, illustrating that the transition is caused by an 
avoiding level crossing and that large lattices sizes are 
necessary to quantify the behavior of this model. The 
first-order nature is apparent in double-peaked magne- 
tization histograms, that were observed for data in the 
"discontinuity" regions for L — 16, indicating a phase 
coexistence. For K/J > 16.3, the spin spin structure 
factor (Eq. Ql) develops Bragg peaks at (7r,7r) (not il- 
lustrated), indicating Neel order. It is interesting to note 
that a similar phase transition between a superfluid and a 
(tt, tt) staggered solid is found in hardcore boson Hubbard 
models with nearest and next-nearest neighbor repulsion 
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FIG. 3: Magnetization (m) and spin stiffness (ps) of the 
superfluid-Neel transition, along cut (B) in the phase diagram 
Fig.Iil Model parameters are h/J = 8, f3J = 320 for L = 16 
and f3J — 400 for the larger lattice sizes. The inset shows 
a double-peaked magnetization probability histogram P{m) 
representing 5 x 10^ Monte Carlo steps at a point on the 
L = 16 data curve in the transition region. The lower peak 
is not Gaussian in shape, as the system is attracted to zero 
magnetization (the half-filled) state. 



Finally we examine the XY superfluid- VBS transition 
along cut (C). As illustrated in Fig. 01 simulation data 



for system size L = 24 does not display an obvious sharp 
discontinuity as in Fig. |3| However, the presence of a 
small discontinuity in m and for L = 32 to 48 is sug- 
gested by the data. The inset of Fig. ^displays a double- 
peaked magnetization probability histogram in the tran- 
sition region, which indicates the presence of a first-order 
phase coexistence. This clearly precludes the existence of 
a continuous quantum phase transition, at least for the 
field value h/J — A that was studied in cut (C) (see 
Fig. GJ . The most immediate conclusion to draw is that 
the superfluid- VBS phase transition is weakly first order, 
either along its entirety (excluding the h = QCP), or 
up to a tricritical point at a field < h < 4. In this 
case, the difficulty in seeing a large discontinuity in the 
superfluid density or plaquette structure factor is due to 
the small h value and the closeness of the magnetization 
to zero. The persistence of a small region of superfluid 
density in apparent coexistence with a finite VBS order 
parameter (for example, the two data points for L = 48, 
K/J = 11.60 and 11.65 in Fig. ^ is due to the first- 
order metastability between the superfluid phase and the 
VBS phase that is obscured by statistical averaging. As 
a check, we observed the Monte Carlo time correlation 
between Sp{TT,0), Sp{0,Tr), and the superfluid density ps 
in the x and y-directions at these points. In fact, we find 
that Ps show no preference for the x or y directions when 
stripe order is present. Rather, both p^ and show very 
strong anti-correlations whenever (Mp)'^ develops Bragg 
peaks. 

In summary, using SSE QMC techniques, we have de- 
termined the ground state phase diagram (Fig. ^ of the 
easy-plane J-K-h model. In addition to the XY super- 
fluid, VBS, and Neel ordered phases observed for = 
13, we observe a large region of fully polarized order, 
which dominates the phase diagram for large h. The 
phase transition to the polarized state is continuous at 
small K/J and strongly first order for large K/J, suggest- 
ing the existence of a tricritical point somewhere on the 
phase boundary for intermediate K/J. Two other phase 
transitions were studied in detail, the superfluid-Neel and 
superfluid- VBS transitions. Both were flrst order for the 
parameter values investigated in detail here. 

As indicated by our data, the J-K-h model does not 
appear to support a region of superfluid- VBS coexistence 
(i.e. a supersolid), which is observed near a similar tran- 
sition between a superfluid and (tt, 0) striped solid phase 
in a hardcore boson Hubbard model 0, Q]. No ad- 
ditional ordered phases were observed in this model, in 
particular, incommensurate VBS stripes (or striped or- 
der away from half filling) which would have been indi- 
cated by Bragg peaks in the q-dependent structure factor 
SpiQx, qy) away from (tt, 0). 

In the context of the h = Q superfluid- VBS transi- 
tion at T = 0, 13 J the existence of a continuous QCP 
does not require a continuous phase transition to develop 
smoothly as h is increased from zero. Conversely, the ex- 
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FIG. 4: Details of the magnetization, spin stiffness and the 
VBS structure factor of the superfluid-VBS transition, along 
cut (C) in the phase diagram Fig. Model parameters are 
h/J = A, I3J = 320 for L = 24 and I3J = 400 for the larger 
lattice sizes. The inset shows a double-peaked magnetization 
probability histogram P{m) representing 3.5 x 10^ Monte 
Carlo steps at a point on the L = 48 data curve in the tran- 
sition region. 



istence of a true continuous phase transition would pro- 
vide additional supporting; evidence for the existence of 
the /i = QCP [121, as well as a further region in which 
to explore the nature of the critical behavior associated 
with the transition from XY superfluid to VBS order. Ul- 
timately, one would like to determine whether this QCP 
is an example of the "deconfined" quantum criticality re- 
cently discussed by Senthil et al. p^ . 

Finally, the inability of any significant region of a spin- 
liquid phase to develop in the vicinity of the QCP mo- 
tivates further searches on related models. Of particular 
interest is the square lattice J-K ring model in a stag- 
gered magnetic field, which could conceivable destabilize 
the superfluid or VBS order near the QCP and promote 
the development of an extended region of disorder. Work 
on this model is in progress ^[ . 
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